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' Abstract. We perform rigorously the charge renormahzation of the so-called 

i-H ' reduced Bogoliubov-Dirac-Fock (rBDF) model. This nonlinear theory, based 

on the Dirac operator, describes atoms and molecules while taking into ac- 
, count vacuum polarization effects. We consider the total physical density 

_ ' including both the external density of a nucleus and the self-consistent polar- 

^ ' ization of the Dirac sea, but no 'real' electron. We show that Ppjj admits an 

asymptotic expansion to any order in powers of the physical coupling constant 
Qph, provided that the ultraviolet cut-off behaves as A --^ e'^'^'^~-^3)/2aph ^ i 
' • The renormahzation parameter < Z3 < 1 is defined by Z3 = Oph/ct where 

- ' a is the bare coupling constant. The coefficients of the expansion of Pph are 

independent of Z3, as expected. The first order term gives rise to the well- 
r~| , known Uchling potential, whereas the higher order terms satisfy an explicit 

recursion relation. 
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1. Introduction and main result 

Renormahzation is an essential tool in Quantum Electrodynamics (QED) [HI [2j 
119] . The purpose of this paper is to perform rigorously the charge renormahzation 
of a nonlinear approximation of QED, the reduced Bogoliubov-Dirac-Fock (rBDF) 
theory that was studied before in [Ml [15l [HI [ITl [HI [12] . This model, based on the 
Dirac operator, describes atoms and molecules while taking into account vacuum 
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polarization effects. It does not need any mass renormalization, hence it is a theory 
simple enough for an investigation of charge renormalization in full detail. 

Before turning to our specific Dirac model, let us quickly recall the spirit of 
renormalization. A physical theory usually aims at predicting physical observables 
in terms of the parameters in the model. Sometimes, interesting quantities are 
divergent and it is necessary to introduce cut-offs. For electrons the parameters are 
their mass m and their charge e (or rather the coupling constant a = e^). Predicted 
physical quantities are then functions F{m,a, A) where A is the regularization 
parameter. Mass and charge are also physical observables and renormalization 
occurs when their values predicted by the theory are different from their 'bare' 
values: 

(1.1) mph = TOph(r7i, a, A) 7^ m and/or aph = aph{m, a, A) ^ a. 

In this case the parameters m and a are not observable in contrast with mph — 
"T'ph("i, a. A) and aph = apii{m, a, A) which have to be set equal to their experi- 
mental values. The relation (|1.1[) has to be inverted, in order to express the bare 
parameters in terms of the physical ones: 

(1.2) m = m(r7iph, aph, A) a = a(r7iph, Oph, A). 

This allows to express any observable quantity as a function F of the physical 
parameters and the cut-off A: 

(1.3) -F(mph, aph. A) = F(m(TOph, aph, A) , a(mph, aph, A) , A) 

A possible definition of renormalizability is that all such observable quantities have 
a limit when A — >■ cxj, for fixed mph and aph- 

Important difficulties can be encountered when trying to complete this program: 

• The physical quantities mph and ofph might be nonexplicit functions of a 
and TO. The corresponding formulas can then only be inverted perturba- 
tively to any order (usually in a) . This is the case in QED 8l [2j [19] . In the 
model studied in this paper we have TOph = m and aph 7^ a, hence only the 
charge has to be rcnormalized. Furthermore aph is an explicit function of 
TO, a and A (see (jl.9p later). Renormalizing our model is therefore a much 
easier task than in full QED. 

• Even when the bare parameters are explicit functions of the physical ones, 
these relations can make it impossible to take the limit A — >■ 00 while keeping 
TOph and aph fixed. As we will explain, in our model (2/37r)aph log A ^ 1. 
To deal with this problem, we let A depend on aph and we investigate the 
asymptotics in the limit aph — 0. 

We now turn to the description of our model. The Bogoliubov-Dirac-Fock theory 
is the Hartree-Fock approximation of QED when photons are neglected [l8l|T7]. The 
associated reduced theory is obtained by further neglecting the so-called exchange 
term. In both models, the system is described by a Hartree-Fock (quasi-free) state 
in Fock space, which is completely characterized by its one-body density matrix P 
(an orthogonal projector for pure states), acting on the one-body space. The state 
P contains both the 'real' electrons of the system (that of an atom for instance) 
and the 'virtual' electrons of the Dirac sea, which all interact with each other 
self-consistently. Therefore, there are always infinitely many particles and P is 
infinite-rank. 
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When the exchange term is neglected, a ground state at zero temperature is 
(formally) a solution of the following self-consistent equation: 



(1.4) 



Here D° = a ■ (— iV) + /3 is the free Dirac operator |27j acting on the Hilbert space 
:= L'^(m.'^,C'^)- For the sake of simplicity we have chosen units in which the 
speed of light is c = 1 and, as the model does not need any mass renormalization, 
we have taken m — 1 for the mass of the electrons. The second term in the formula 
of D is the Coulomb potential induced by both a fixed external density of charge v 
(modelling for instance a smeared nucleus) and the self-consistent density 
of the system (see below). In (|1.4p . a is the bare coupling constant that will be 
renormalized later and /iG (—1,1) is a chemical potential which is chosen to fix the 
desired total charge of the system. We have added in (|1.4|) the possibility of having 
a density matrix 0^5^ X{/i}(-^) the Fermi level, as is usually done in reduced 
Hartree-Fock theory [26]. So the operator P is not necessarily a projector but we 
still use the letter P for convenience. Later we will restrict ourselves to the case of 
P being an orthogonal projector. 

Equation (|1.4p is well-known in the physical literature. A model of the same 
form (including an exchange term) was proposed by Chaix and Iracane in p]. Also, 
similar equations are found in relativistic Density Functional Theory, usually with 
additional empirical exchange-correlation terms and classical terms accounting for 
the interactions with photons, see, e.g., [TTl Eq. (6.2)] and (THl Eq. (62)]. Dirac 
already considered in [7] the first order term obtained from (jl.4l) in an expansion 
in powers of a. 

Let us now explain the exact meaning of pp -1/2- The charge density of an 
operator A : ^ Sj with integral kernel A(x, y)a,a' is formally defined as pa{x) — 
Ylt=i ^(-^i — Trc4(A(a;, x)). In usual Hartree-Fock theory, the charge density 
is pp{x). However, as there are infinitely many particles, this does not make sense 
here. In (|1.4p . the subtraction of half the identity is a convenient way to give a 
meaning to the density, independently of any reference. One has formally 

Pp-i/2{x) = P p-p± jx) = |<^r(2;)P - \'Pt{x)\^ 

where {<Pj~}«>i orthonormal basis of Pio and {'P^Ji^i is an orthonormal basis 
of (1 — P)Sj. As was explained in [TB], subtracting 1/2 to the density matrix P of 
the Hartree-Fock state makes the model invariant under charge conjugation. 

When there is no external field, 1^ = 0, Equation (|1.4p has an obvious solution for 
any /i G (— 1, 1), the Hartree-Fock state made of all electrons with negative energy: 

P = P° :=X(-oo,o)(^°), 

in accordance with Dirac's ideas [11[5J|B]. Indeed Ppo ^1/2 = 0, as is seen by writing 
in the Fourier representation 

a. ■ p + P 



{P^_ - l/2){p) 



and since the Dirac matrices are trace-less. This shows the usefulness of the sub- 
traction of half the identity to P, since the free vacuum now has a vanishing 
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(1.6) 



density. For a general state P, we can use this to write (formally): 

(1.5) PP-1/2 = PP-1/2 - PpO-l/2 = Pp-p° ■ 

When P belongs to a suitable class of perturbations of P° (for instance when P—P^ 
is locally trace-class), the density pp_po is a well-defined mathematical object. We 
will give below natural conditions which garantee that P — P° has a well-defined 
density in our context. 

In the presence of an external field, ly =/= 0, Equation (jl.4l) has no solution in 
any 'reasonable' Banach space |15| and it is necessary to introduce an ultraviolet 
regularization parameter A. The simplest method (although probably not optimal 
regarding regularity issues [12]) is to impose a cut-off at the level of the Hilbert 
space, that is to replace i5 by 

{/ e L\M.^;C^), supp(/) C P(0;A)} 

and to solve, instead of ()1.4p . the regularized equation in 

j D = TIk[d° + a{pp_po - 1^) * Ha 

where IIa is the orthogonal projector onto Sja in = L^(R'^;C"'). 

Existence of solutions to (|1.6p was proved in [T3] for fi — and in "TT for /i 7^ 0. 
The precise statement is the following^: 

Theorem 1 (Existence of self-consistent solutions to (|1.6p . [Ullll]). Assume that 
a ^ 0, A > and /i G [—1, 1] are given. Let v in the so-called Coulomb space; 

C:= (/ : / \k\-^\m\''dk<^ 

Then, Equation (jl.6|) has at least one solution P such that 

(1.7) P-P^ e ©2(i3A), P£(P - P")P£ e 6i(i3A), pp_poeCnL\M?). 

All such solutions share the same density pp_po . 

In (II. 7p . &i{Sja) and &2{^a) are respectively the spaces of trace-class and 
Hilbert-Schmidt operators [25 on and P° = 1 — P°. The method used in 
[151 [T^j was to identify solutions of (|1.6p with minimizers of the so-called reduced 
Bogoliubov-Dirac-Fock energy which is nothing but the formal difference between 
the reduced Hartree-Fock energy of P and that of the reference state P". Note 
that due to the uniqueness of pp-po the mean-field operator D is also unique and 
only S can differ between two solutions of (|1.6p . 

Let us mention that it is natural to look for a solution of (|1.6p such that P — P" 
is a Hilbert-Schmidt operator on Sja- If P is a projector, the Shale-Stinespring 
theorem [SJ then tells us that P yields a Fock representation equivalent to that of 
P^. Even when P is not a projector, it will be associated with a unique Bogoli- 
ubov mixed state in the Fock space representation of P^. This is a mathematical 



^To be more precise, in |12l Theorem 1] , only the existence and uniqueness of minimizers of the 
reduced BDF functional are stated. Elementary arguments based on convexity allow to deduce 
Theorem [1] from the results of |12| . 
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formulation of the statement that P should not be too far from P° . Indeed, if P is 
an orthogonal projector, one has (see [111 Lemma 2] and [Ml Lemma 1]) 

p-p"_e I ^ po (p _ pO)po ^ ^^^^ ^ c n L\m'), 



therefore, in this case, (jl.7l) is just equivalent to the Shale-Stinespring condition 
P-P° e 62(-$5a). 

The property P^{P - P°)P^ e 6i(^a) allows us to define the total 'charge' of 
the system by (see [H]) 

Trpo (P - P°) Tr P°(P - P°)P° + Tr P°(P - P!?)P°. 

When P is a projector, the above quantity is always an integer which is indeed 
nothing but the relative index of the pair (P, P°), see [Mid]. Varying ^ allows to 
pick the desired total charge. Indeed, if v is small enough and /i = 0, then one has 
II -P — P- II < 1 and the relative index vanishes: Tr^o (P — P") = 0. 



It is very important to realize that solutions of (|1.6p are singular mathematical 
objects. This fact is precisely at the origin of charge renormalization. In [T2l 
Theorem 1], the following was proved: 

Theorem 2 (Nonperturbative charge renormalization formula [12]). Assume that 
a ^ 0, A > and ^ € (-1,1) are given. If ly € C n L'^(R^), then pp_po G L1(R3) 
and it holds 



(1.8) Pp-P[ 



/ i/-Trpo(P-P°) 



1 + aPA 

In this formula. Pa is an explicit function of the ultraviolet cut-off A (see the 
comments after ()2.4p and ()B.10|) ). which behaves like 

i^A = flogA-f + ^ + 0(l/A^). 

dTT 97r JTT 

Let us emphasize that (11.81) is non perturbative and holds for all a 5^ and all 
^ E (—1, 1). Theorem [5] shows that the operator P—P^ is in general not trace-class: 
if p _ po g 6i(i0A), then it must hold Tr^o (P - po ) = Tr(P - P°) ^ J^,, Pp-po . 

In our model we have two possible definitions of the charge of the system: jj^g v — 
Trpo (P — P°) and ji^^ii^ — pp-po )■ In practice it is the electrostatic field induced 
by the nucleus (together with the vacuum polarization density) which is measured, 
hence it is more natural to define the charge by means of the density. By (|1.8p . the 
total Coulomb potential is, at infinity, 

, 1 /r3 (i^ - Pp-p« ) T+fsI ( /r3 ^ " Trpo (P - P"_)] 
a(y — pp_pa ) * - — - ^ a j — j — ^ ■ — ■ —. 



Let us assume for simplicity that we put in the vacuum (/i = 0) a nucleus containing 
Jj^3 V = Z protons and which is small enough in the sense that <C 1. Then 
Trpo (P — P") = by [TS] Theorem 3] and we see that at infinity the potential 
induced by the nucleus is not aZ /\x\ as expected, but rather a-p\iZ /\x\ where 

(1.9) ttph = ^3" , with Z3 = , ^ „ = 1 - aph^A- 

i + Qf-DA 
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The charge renormaUzation constant Z3 is weh known in QEE@[8l[2l[l9]. The value 
of a is not observable, aph is the real physical constant since we always observe the 
nucleus together with the vacuum polarization density. Its experimental value is 
Qfph 1/137 . 

In our theory we must fix aph and not a. Using (jl.9p we can express any physical 
quantity in terms of aph and A only. Unfortunately it holds aph^A < 1 hence it 
makes no sense to take A — > c» while keeping aph fixed (this is the so-called Landau 
pole [5D]) and one has to look for a weaker definition of renormalizability. The cut- 
off A which was first introduced as a mathematical trick to regularize the model 
has actually a physical meaning. Because of the above constraint aphi?A < 1, a 
natural scale occurs beyond which the model does not make sense. Fortunately, 
this scale is of the order e'^'^/^"p'>, a huge number for aph — 1/137. 

It is more convenient to change variables and take as new parameters aph and 
Z3 = 1 — aph-BA, with the additional constraint that < Z3 < 1. The new 
parameter Z3 is now independent of aph and the natural question arises whether 
predicted physical quantities will depend very much on the chosen value of < 
Z3 < 1. The purpose of this paper is to prove that the asymptotics of any physical 
quantity in the regime aph <C 1 is actually independent of Z3 to any order in aph, 
which is what we call asymptotic renormalizability. Note that fixing Z3 £ (0, 1) 
amounts to take A ~ (7e3'^(i-^3)/2Qph i 

Instead of looking at all possible physical observables, it is convenient to define 
a renormalized density pph. Following [15! , we define it by the relation 

(1.10) aphPph ^ a{v - pp_po) 

in such a way that D = — aph/Cph * Nl"^- This procedure is similar to wave- 
function renormalization. By uniqueness of pp^po we can see pph as a function of 
ttph, V, p and A (or Z3). For the sake of clarity we will not emphasize the depen- 
dence in v and p which will be fixed quantities. Also we will use the same notation 
Pphictph, A) or Pphiciph, Z^), depending on the context. The self-consistent equation 
for Pph was derived in |15) and it is mentioned below in Section [2] 
From now on, we will assume that 

p^O. 

For small external densities v, this means that we will be looking at the vacuum 
polarization in the presence of the nucleus, without considering any real electron 
(that is, pph is the renormalized density of the nucleus containing both the bare 
density v and the vacuum polarization density pp_po). We will explain in Section 
[2] that one can expand pph = /Oph(Q^ph, A) as follows: 

00 

(1-11) Pph(aph, A) = y^(aph)"t^n,A 

n=0 

where {i/„,a}„ C L'^{R^)nC is a sequence depending only on the external density 
ly and the cut-off A. This sequence is defined below in Section [2] The series (|1.11|) 
has a positive radius of convergence, which is however believed to shrink to zero 
when A ^ 00. 

Assuming 9 decays fast enough (see condition (|1.12p ). we will prove that for any 
fixed n, the limit Vn.A i^n exists in L'^(U.^) n C. This is what is usually meant 



'The renormalization constant Z3 should not be confused with the nuclear charge Z = u. 
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by renormalizability in QED: each term of the perturbation series in powers of the 
physical ctph has a hmit when the cut-off is removed. The sequence {t'njn is the 
one which is calculated in practice [H [T31 El [TO] ■ One has for instance 1^0 = 1^ and 



vi * \x\ 



1 

3^ 



di(t2- 1)1/2 



2 

i2 



1 



-2\x — y\t 



x-y\ 



dy, 



the Uehling potential \28\ 123] . All the others can be calculated by induction in 
terms of vq, Vn-i, as is explained below in Section [2j 

The next natural question is to understand the link between the well-defined, 
cut-off dependent, series (jl.lip and the formal series X]^S^o('^ph)"^n- Recah that 
<XphB\ < 1 by construction, so it is in principle not allowed to take the limit A — > 00 
while keeping aph fixed: we rather want to think of Z3 = 1 — aphi?A as being fixed. 
The main result in this paper is the following 

Theorem 3 (Asymptotic rcnormalization of the nuclear charge density) . Consider 
a function v G L^(M'') n C such that 



(1.12) 



log(l + |fc|)2^+2|P(A:)|^rffc < 00 



for some integer N. Let Pph(Q^ph7 -^3) be the unique physical density defined by 
(jl.lOP with fi ^ 0, = 1 — Q!phi?A O'lT'd Qfph = Z^a. 

Then, for every < e < 1, there exist two constants C{N^e,v) and a{N,e,i'), 
depending only on N , e and v , such that one has 



(1.13) 





N 






Pph(aph, Zz) - ^ t'„(aph)" 


s=:C(7V,e,^) a^+i 




n=0 


L2(R3)nC 



for all ^ aph ^ a(A^, ^^ v) and all e ^ Z3 ^ 1 — e. 

The interpretation of Theorem [3] is that the renornialized density Pph(Q^phj ■^a) 
is asymptotically (meaning up to any fixed order N) given by the formal series 
Sn>o(*-*^ph)"^ni uniformly in the rcnormalization parameter Z3 in the range e ^ 
Z3 ^ 1 — e. Therefore, for a very large range of cut-offs, essentially 

the result is independent of A and it is given by the formal series X]n>o('-'^ph)"^n- 
Our formulation of renormalizability is more precise than the requirement that each 
t^n,A converges. It also leads to the formal perturbation series in a very natural way. 

A natural question is to ask for the convergence of the perturbation series 
Sn>o('^ph)"^n- It argued by Dyson in [5] that it is probably divergent, but we 
are unable to transform his argument into a rigorous mathematical proof. We will 
make more comments on the series X]n>o('^ph)"^n ^^'^ next section. 

Remark 1. We will provide explicit formulas for the sequence {vn} later in Section 
We will in particular see in the proof that under Assumption (|1.12p . one has 
Vn G i^(R"^) n C for all ^ n ^ N . Therefore the approximation series of order N 
appearing in (|1.13p . X)^=o('-'^ph)"^nj '■^ well-defined function of L'^(R^) nC. 

Remark 2. The space L'^iM.^) CiC is the natural space which occurs in this theory. 
In particular the Coulomb norm is nothing but the classical electrostatic energy 
which appears in the reduced BDF energy functional. Our result can be extended to 
Sobolev spaces H^(U.^) provided v is smooth enough. 
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Remark 3. It would be interesting to extend this result to the case of atoms with 
'real' electrons. This amounts to taking /i sufficiently close to 1 at the same time 
as ckph is small. However this case is more difficult than what is done here: an 
additional expansion of the electronic charge density in powers of aph is needed. 

The proof of Theorem [3] (given m Section [5] below) is divided into two steps. We 
first estimate the difference (see Lemma [5. ip 



(1.14) 



Pph(aph, A) 



N 

E 

n=0 



^'n,A(aph)'^ 



s$ Ci(iV,e,i^)(aph) 



N+l 



L2(R3)nc 



for a constant Ci{N, e, v) depending only on N, e and i^, and under the assumption 
that e ^ Z3 = 1 — aphBA ^ 1 — e. This amounts to expanding the solution of the 
self-consistent equation (II. 6|) up to the iVth order in aph while controlling the error 
term uniformly in A. Then we show in Lemma l5.2l that 

II II <r ^^(iV,^^) 

ll^'n.A |^2/"<<N"'^ ^ 



(1.15) 



VO s$ n s$ iV, 



for a constant C2{N, v) depending only on N and v, leading to the bound 



(1.16) 



N 



i^n,A(aph)" 



N 

E 

n=0 



< ("ph 



e^+i(l-e) 

since by assumption (_Ba)~^~^^" ^ («ph/e)^''~^~"- The main result then follows 
from (|1.14p and p. 161) . All these bounds strongly use the explicit recursion relations 
defining the sequences {vn,A} and {v,i\, as well as tedious estimates on the nonlinear 
terms appearing in these relations. 

The rest of the paper is organized as follows. In Section[2]we define the sequences 
{i'n,K\ and {vn\ by their respective recursion formulas and we discuss some prop- 
erties of the latter. In particular, in Theorem 21 we give a simple estimate on 
ll^"llL2(jj3)nc- In Section 13] we present estimates on the different terms appearing 
in the recursion formulas. Of particular interest will be the density i^i.a giving rise 
to the Uehling potential. In Fourier space, we have ?i,A(fc) = U\{k)v{k) for an 
explicit function UA{k) which is studied in Section [ST] The proofs of Theorems H] 
and [3] are respectively provided in Sections HI and [SJ Some other technical proofs 
are provided in Appendices \^ [B] and [C] 

Acknowledgment. The authors are grateful to Christian Brouder for interesting 
comments. M.L. would like to thank Jan Derezinski and Jan Philip Solovej for 
stimulating discussions. 

2. The two sequences {vu.a} i^n} 

In this section we derive formulas for {i'n,A} and {vn}, and wc make some com- 
ments on the latter. 

2.1. Definition of {fn,A} and We start with the self-consistent equation 

()1.6p with cut-off, assuming /i = 0. Note that in the regime of interest in Theorem 
|3l we have a = aph/Z^ ^ ctph/e. When < 1, it is known that 

^ o-{D) hence ^ = in (jl.6|) . see [EJ Theorem 3] and Lemma 11]. Therefore 
assuming a(i/, iV, e) ^ e{TT^/^2^^/'^ g) "'^ in Theorem[3l we automatically have that 
6 = and F = is unique. 
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The idea is then to expand the self-consistent equation 

(2.1) P - X(-oo,o)nA(i?° + «(pp-po - ly) * \x\-')Ua 

in powers of a by means of the resolvent formula. This method was already used 
in [H] to prove existence and uniqueness of solutions. We define 



(2.2) F„,a(mi: •■•7/^fc) P 



1 /■°° 1 -A- / 1 1 



drj 



where we recall that IIa is the orthogonal projector onto i^A in L'^(M'^,C'') and 
^i, ...j^n G C. We will always use the simplified notation F„_a(/i) :— Fn,K{p-, •■•,m) 
and Vf^ := J-^^{i''^b{q.2A))- Note that by Furry's theorem F2j,A = for all j, see 
[m p. 547]. We also introduce 



F^{^J) :=^F„,a(^). 



n^3 



The self-consistent equation (12. ip may then be written in terms of the density 
in Fourier space [HI [T5] , as 

(2.3) p]r^o{k) = ~aBAik){p^r~Po{k) -VA^k)) + FA{a{iy - pp_po)) 

where the function BA{k) is given by 

27: Jo VTTA^-\k\z/2 

with ZA{r) = + A2- ^1 + (A - r)2) /r, see [12]. The formula for Ba(/c) is well- 
known (but in most previous works the second term was ignored, see for instance 
[21]). 

Defining C/a(|A:|) ^ Ba- BA{k) where Ba = Ba(0) and UA{\k\) Ba with 
C/a(2A) — _Ba, we get the renormalizcd equation 

(2.5) (l - aphC/A)pph + ^A(aphPph) = va 

with the renormalizcd coupling constant ofph ct/(l + olBa) and the renormalizcd 
density aphPph = ci{u — pq) (see [15]). For convenience, we will denote by IAa the 
operator of multiplication by the function [/a (1^1) in the Fourier domain. Hence we 
can write the self-consistent equation (|2.5p in direct space as 



(2.6) (l - aphWA)pph + ^A(aphPph) = I'A- 



We now expand the unique solution pph = Pph(aph, A) of (j2.6p in powers of ttph- 
Writing a formal series 

(2.7) = ^(aph)"i^„,A 
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we find that the functions Vn^K must satisfy the following recurrence relation 
t'o.A = vk. 



(2.8) <^ 



Vn,K=UKVn-l,h+'^ ^ -Fj, A (i^ni , A, ■ • ■ , J^n^ ,a) , Vn ^ 2. 

i=3 TiiH \-nj — n-j 



Note that the operator Ufi^ is bounded by U (2A) = i?A on and that, as we will 
see later in Corollary 13. 1[ each F^- a is continuous on with values in nC The 
sequence {m„.a} is thus well-defined in nC. Using estimates from [T3] it can be 
proven that the series (|2.7p has a finite radius of convergence in nC, but this is 
not needed for the moment and we can stay at a formal level in this section. 

We can now formally pass to the limit as A — )• oo and define by induction a 
sequence {i^n} by 



(2.9) 



^0 



vi =Uv, 



j—3 niH \-nj—n—j 



where the Fj are defined similarly as the Fj^\ with IIa removed and U is the operator 
of multiplication by the function t/(|fc|) in the Fourier domain, defined by 



(2.10) ^M-ii^mf/A(.) = ^/^^^^ 



dz 



12 - 5r2 



[r^ - 2) log 



V4~- 



v4 + r^ — r 



2.2. On the series {vn}- The recursion formula (j2.9p defining {vn} contains two 
terms. The first term Uvn-i is a simple multiplication operator in Fourier space, 
by the function C/(|fc|) which diverges at infinity. The second term involves the 
nonlinear functions Fj's. If only the first term with U were present, the series 
Sn>o '^"(o^ph)" would only converge when the Fourier transform v has a compact 
support, the radius of convergence depending on the size of this support. If only the 
nonlinear terms were present, the series would have a finite radius of convergence 
by the estimates of T3] and of Section [X^ 

However when the two terms are combined, the situation is much more compli- 
cated. The nonlinear terms act like convolutions in Fourier space, hence even if P 
has a compact support in the Fourier domain, the support of ;?„ will probably grow 
with n. A careful study of the mixed effect of the multiplication by the divergent 
function U and the nonlinearities seems rather difficult. We will prove the following 
estimate: 
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Theorem 4 (Estimate on {vn}n^i)- There exist universal constants A and K such 
that 

(2.11) ||(l+Z^)'"-"^„|L.nc 

^ A"+imax|||(l+Z^)"V|U.nc , {Klog{m))'^\\{l+Ur'i^\\lt^c 

for all A ^ 1, m Cz N and ^ n ^ m. 

Even if we assume that v decays fast enough in Fourier space, for instance 

Vn^O, ||(l+Z^)V||^.nc^C", 

the above estimate (|2.1ip does not imply that the series X]n>o '^"('^ph)" 
vergent for aph small enough. Although our estimate (j2.1ip is certainly far from 
optimal, as we have already mentionned, it is expected that the series does not 
converge in any appropriate sense [9]. 

It is sometimes argued that the series could be Borel summable. The Borel 
transform is defined by 

Bit) = y: 

If B(t) is a convergent series (for an appropriate norm) having a holomorphic ex- 
tension to a domain containing the positive real line, such that 

B(aph) := / B{t)e "p" dt 
Jo 

makes sense in an appropriate neighborhood of ctph = 0, one may see B{aph) 
as the physical density, whose series J2n^o('^pt^)"^'^n is only asymptotic. Proving 
such results mathematically is hard, even for the model studied in this paper. Our 
estimate (|2.1ip does not even allow to define the Borel transform B{t) in L^(R'^)nC. 

But Borel summability is not the only tool to construct a physical density pro- 
viding the correct asymptotic series. For the model studied in the present paper, 
we have several natural families of functions of ckph, the cut-off densities 

(2.12) Pph(aph,Ce3(^-^^)^/'"^'0 

obtained by minimizing the reduced BDF energy with a cut-off A = (763(1-^3)^/20^^ 
and using the relation (jl.lOp . Each such density (|2.12p has (for fixed C and < 
Z3 < 1) the required asymptotic series in aph by Theorem |3l and it solves the 
self-consistent equation (|1.6p with the corresponding cut-off A. Furthermore this 
solution has the benefit of being well-defined even when aph is not small, allowing 
for the description of nonperturbative physical events. 

The rest of the paper is devoted to the proofs of Theorems [3] and S) 

3. Some preliminary results 

In this section we state two preliminary results that will be useful in the proof of 
our main results. Theorems [3] and SI The corresponding lengthy calculations will 
be provided later in Appendices El [B] and [C] 

Notation. In the whole paper we use the notation E{r) = (1 + |rp)i/^ for r e 
or r e M. 
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3.1. The Uehling multiplier U. The operator Z^, defined previously as the mul- 
tiplication by the function U in the Fourier domain, plays a major role in the 
definition of the sequence {i^n}- In this section, we provide precise estimates quan- 
tifying the convergence of U\ towards U when A — >■ cx), which will be very useful 
in the proof of Theorem [3l 



Proposition 3.1. Let A ^ 1 and denote by Ua, the function defined on Rj^ by 
(3.1) C/A(r) = 



BA-SA(r) when ^2A, 
otherwise. 

Then, for all r G it holds liniA-i-oo C^a(^) = U{r). Moreover, for kq ~ 157r/2, 



(3.2) Vm ^ 0, 



Ua-U 



m+l 



max 



{(l + B 



Finally, one has for a universal constant ki (given in Lemma \B.l\ below) 
(3.3) V0==:r<2A, s$ C/A(r) sC ki (l + t/(r)) . 

Proposition 13.11 is proved in Appendix |B1 Note that the uniform estimate 
will later yield our estimate (ll.lSp on i^n.A — i^n (see Lemma [5?2|) . More properties 
of U and Ua are provided later in Appendix |^ 

3.2. The nonlinear terms Fn,A and F„. In this section, we provide estimates 
on the functions Fn.A and Fn, which will be one main ingredient in the proof of 
Theorem [31 We recall that F2n,A = F2n = by Furry's theorem (see [HI p. 547]). 
In order to state our main result, we introduce the functions 



(3.4) i^^,A(M)=P 

l^T^ J-oo ^ "I 

for any ti ^ 3, = (^i, • • • ,^„) E C" and e = (ei, • • • ,e„+i) E {-1,0, 1}"+^ Here, 
we have used the notation 

(3.5) n^'^ Ha, H^"'' := 1 - Ha and nj^°^ := 1 = H^'^ + u'^'l 
The main result of this section is the following 

Proposition 3.2 (Estimates on F^ a)- Let m E N, A ^ 1 and e E {-1,0, 1}"+^ 
Assume that n ^ 3. Then, there exist universal constants C and K such that 

(3.6) 1(1 + w)"^;u(M)IL.nc ^ In J J n 11(1 +")"M.llc > 

J=l 

for all /i = (/ii, • • • , /i„) E C". Here, n{e) — 1, if at least one ej is equal to —1, and 
n{e) = otherwise. 

By ([22]), ^ and we can write i^i^^ '^^ = K,A and = K- 

Therefore the following is a byproduct of (13.61) : 

Corollary 3.1. Let m ^ 0, A ^ 1 and n ^ 3 an odd integer. Then, 

(3.7) max{ 1(1 +Z^)™F„,A(M)lli.nc ' 1(1 + U)"^ F^l^.^^ } 

n 

«;C"(iflogn)™ni(l+^)'>jlc' 
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for any /i — (/ii,--- , fin) G C". Here, C and K refer to the universal constants 
given by Proposition \3.2[ In particular, the functions Fn,A and _F„ are continuous 
on C" with values in L'^ f] C. 

Recall — F2k.K = hence only the case of n being an odd integer is relevant. 
The estimates of Proposition are an adaptation of ideas of (14) . in which similar 
bounds were computed (see, e.g., Lemmas 15 and 16 in [2]). Notice however that 
the projector IIa was never mentionned in |14j since A was a fixed number. We 
focus here on the limit A — > cxd and we need to quantify the dependence on A of 
the estimates on the functions ^. The proof of Proposition l3.2l is provided below 
in Appendix [Cj The factor {KXogn)'^ comes from (jA.7|) of Lemma [A. 31 and the 
constant K is also the one appearing in Theorem |4l 

4. Proof of Theorem [4] 

This section is devoted to the proof of our estimate (|2.11[) on the nth order 
density Vn- The definition of Vn^A being very similar to that of Vn, our proof also 
provides the following 

Proposition 4.1 (Estimates on Vn^A)- There exists A> Q such that 

(4.1) \\{i+ur-un4^,^^ 



^ A 



n+l 



max 



{ II (1 + z^rnu^nc, (i^ log(m)) II (1 + z^rni } 



for any A ^ 1, m G N and ^ n ^ m. 

We postpone the proof of Proposition 14.11 and first complete that of Theorem U) 

Proof of Theorem^ We split the proof into three steps. First, we estimate by 
means of (|3.7p . the following norms: J,„^„ := ||(1 + U)"^^^'' I'nW L^nc ■ 



Step 1. Let m e N and denote 



n=0 



The polynomial P„i (t) satisfies for any t ^ 

(4.2) P„(t) {l + t + t')\\{l+Ur,y\\L2nc + QmitP^{t)), 

where (C and K are the constants of Provosition lST^) 



(4.3) QM := u + J2C^K log jY 

3=3 



Let us assume first that n = 0, 1, 2. By (|2.9p . we then have Vn = U^v, hence 
(4.4) Vn = 0,l,2, Jrn,n = \\il+Ur-^U^v\\L-^c^\\(l+Urv\\L^nc- 
We now turn to the case n ^ 3. By (|2.9p . we have 



J2 E ii(i+z^)'"-"f^2,+iK, 



I '^n2j + i)\\L^nC, 



3^2j+l^n 2j + l 

J2 nk=n-2j-l 

k=l 
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hence, by Corollary 13.11 

J™,„ ^ ||(l+Z^)"~"+V„_i|U2nc+ J2 E ^''^'^ 

3<2j"+l^n 2j + l 

J2 Tlfc=Tl-2j-l 

;c = i 
2j+l 

X (xiog(2j + i))™-" n ii(i+z^r-v„jic. 

k=l 

Since ||(1 +i/)"""i/„J|c ||(1 + W)"'""'!^nJ|L2nc, we arrive at the inequality 

i 



(4.5) J„,„ Jra,„-l + C^ (A'logj)™-" ^ ( n J, 



i=3 J ^ fc=i 

k = l 



Combining (|4.4I) with (14.51) . we obtain 

/'m(i) ^ {l+t + t^)\\{l+Uri^\\L2r,c+tP.m{t) 



m n / ^ \ 

J2 E CH^iKlogjY^- n ) • 

n=3 j=3 i ^ k=l ' 

nk=n-j 

k = l 



By Fubini's theorem, it holds 

(4.6) P,n{t) < {1 + t + t^)\\{l +Uy"'iy\\L2nc + tPm{t) 



]=Z p=0 j ^ k=l 

12 nk=p 

k = l 



+ Y cH\K\ogjr-= E E ( n J-^^r.kt'-'' ) ' 

Noticing that 



E E [EJrn.^kt"'') ^Prnity, 



p=0 i 

k = l 

we deduce (|4.2p from (|4.6I) . This completes the proof of Step [TJ In the second step 
of the proof of Theorem |4] we compute suitable bounds on Qm near the origin. 

Step 2. Let m ^ 3. There exists a positive constant A{C,K), depending on C and 
K , hut not on m, such that 

(4.7) Q„,(m) < 2u, for any < w U,n ■= 



{K log(m))™/2' 
By the definition (|4.3p of Q™, we have 

rn rn ^ ■ 

Qm{u) sc: u + 5]C^"(if log(j))"-^«^ <^u+{K\og{m)rJ2{lFr^y' 
— ' ^ — ' \K log TO/ 

J=3 j=3 ^ 

hence when 2Cu ^ iflogm and 2C^(inog(TO))"-3y2 j^^j^g Q™(-u) s$ 2-it. 

This ends the proof of Step [2l 

Step 3. Conclusion of the proof of Theorem^ 
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Since the coefficients Jm,n are non negative, the function t i— > tPm(t) is cither 
identically equal to (then (|2.11|) is straightforward), or increasing on R+. In the 
second case, it tends to oo as t oo, and there exists a unique Tm > such that 

(4-8) ^rn-Pjni'-^m) — ^m- 

Two situations may then occur. If ^ 1/4, by (|4.2|) and (|4.7p . 

Pm{t) < 2\\{l+U)''^i^\\L2nc + 2mnW 
for aU t 1/4. Hence P™(t) 4||(1 +Z^)"i/||i2nc and 

Jm,n ^ 4"P™(i) < 4"+l||(l +W)'"i.|U2nc. 

Otherwise ^1/4 and in this case we can deduce from (|4.2p . (14. 7p and (14. 8p that 
J/m/T™ ^ 2|1(1 +Z^)"j.||L2nc + 2C/,„. This gives 

J- rn ". ^ 



4||(l+Z^)"V|U2nc 
Combining with (|4.8I) again, we are led to 

Tj An+l 

^11(1 +^rHr2t^c- 

Estimate (pTTTj) then follows from (HT7)) . □ 
We now turn to the proof of Proposition |43] 



Proof of Proposition \4-l\ The proof is almost identical. Denoting „ := ||(1 + 
U)"^~'^i'^j^\\]^2^c and introducing the polynomial function P^{t) given by 

m 

we deduce from the definition (|2.8I) . and from (13. 3p and (13. 7p that 

(4.9) P^it) ^ (1 + «it + niem+Uy^^uU^nc + ('^i - l)i^'™(0 + Q,n(tP,^(t)), 

for all i 0. Estimate (|4.ip then follows by applying to (14. 9p the arguments of 
Steps H and [3] of the proof of Theorem [H □ 

5. Proof of Theorem [3] 

This last section is devoted to the proof of our main estimate (|1.13p . The proof 
relies on the identity 
(5.1) 

Pph — VA + "phC^APph - ^ ap5J+-^_F2ji+l,A(Pph, ' ' ' , Pph) — Olp\^^GN+l,A, 
3^2n+1^7V 

where we denote 



(5.2) Gw+i,a:-P 



1 1 



27r 13° - QfphPph * I • I ^ + 17? 

Af + l 



X 

3 



V + 1 \ 

n(nA(Pp.*^)nA^^)d7?j. 
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The formula (|5.ip fohows from Cauchy's formula applied to (|2.6p . As mentioned 
in the introduction, the proof of (|1.13p naturally splits into two steps: we first 
establish that, under the assumptions of Theorem |31 the error term 

N 

(5.3) i?7v(aph, A) := pph(aph,A) - ^ i'„,Aaph- 

n=0 

is controlled by a factor a^^^^ (up to some multiplicative constant depending only 
on N , V and e). In a second step we estimate the differences i'n,K — i^n and deduce 
(|1.13p . More precisely, the remainder satisfies the following 

Lemma 5.1. Let N E N and < e < 1. Assume that e ^ Z3 = 1 — aphi?A ^ 1 — e 
and Af^ ■— \\{^ +i^)^~^^'^\\L^nC < Then, there exist two constants C{m,e,Afj\i) 
and a(iV, e, TV^v); depending only on N, e and Mn, such that 



N+l 
ph ' 



(5.4) i?w(aph,A) ^C{N,e,NN)a 

L^nc 

for all Qfph a(iV, e. A/at). 

As for the differences j^„,a — i^n, we have the 

Lemma 5.2. Let A ^ I and N € N. Assume that Af^ -.^ ||(l+Z^)^+V||i2nc < 00. 
Then, there exists a constant C{N,AfN), depending only on N andNN, such that 

(5.5) Vn,A - l^n ^ 

for aZ/ < n < TV. 



L^nc (1 + Sa)^+i-"' 



Combining Lemmas 15.11 and 15.21 we can complete the proof of Theorem [3l 

Proof of Theorem\^ Our assumption (|1.12p (together with (jA.3[) ) means that Nn '■= 
11(1 +Z^)^+ij/||i2nc < 00. It follows from that 

N N 

Pph(aph, A) - ^ J^„(aph)" = ^Ar(aph, A) + ^ (z^„,a - i'„)(aph)". 

n=0 n=0 

Hence by and (|5?5t . 
(5.6) 

N ^ rv" 

n ^ ,-(/AT Kr \ N+l , ^/AT «r \ ph 



n=0 



Pph(aph) - V vn<^ C(A^, e, A/-w)a;^+i + C{N,Nn) 



n— ^ ' 



for any number ctph sufficiently small. In our setting we have i?A ^ e/o^ph and the 
result follows. □ 



It therefore remains to show Lemmas 15.11 and 15.21 
Proof of Lemma I5.il Let us introduce the notation 
(5.7) rAr(aph) := (aph)"^"^i?Ar(aph). 

We want to establish a bound on rjv independently of aph. By (j5.3p . this requires 
to estimate pph and i^^^A (which was already done in Proposition 14. 1 1) . 

The first step of the proof will be to bound pph independently of aph. Let us 
recall that a ground state for the reduced Bogoliubov-Dirac-Fock model satisfies 

llttphPphllc = \\a{pQ - v)\\c < a\\v\\e 
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(see [m Eq. (33)]). Since ctph = ^30, this provides 
(5.8) \\ppi,\\c Z-'Mc ^ e-'Mc. 

Note that we however do not have any a priori bound in L^(R^). Inserting (|5.3 



and (|5.7p in (jS.ip and using ()2.8p . we get 



(5.9) Tat = ttph^ArAr + Uai^na + Gn+i^a + ^ " " x 



Af(Ar+2) 

ph 
fc=Ar+l 

X ^ ^ -F2n+l,A(Wpi, ■ ■ • ,Wp2„^J, 

3^2n+l^Af piH |-p2n + l=fc-2n-l 

where ujp = Vp,A for ^ p ^ iV, and w^v+i = fw. It rests to estimate all the terms 
of the right-hand side of (|5.9p . 

For the first term, we recall that aphlC^Al ^ Q^ph^A = 1 ^ .^3 ^ 1 ^ e, therefore 

(5.10) llaphZ^ArArllL^nc < (1 - f)\\TN\\L^nc- 

The second term can be controlled by using (|3.3p and (|4.ip . which provide a positive 
constant C(iV,A/jv), depending only on iV and Nni such that 

(5.11) \\UAyNA\L-nC ^ Ki||(l +Z^)l/^v,A||L^nC 5S C(iV, A/^jv). 

As for the function Gjq+i A, we first recall that 
(5.12) 



(1 _ IHIc ) \D^\ < \d' - «phPph * ^1 ^ (1 + ^«Ph IIHIc ) 

for all aph < 7r-i/f^2-"/f^e ||;/||^^ (see [HI p. 4495]). Hence, the operator D° - 
ttpiiPph * I • 1^^ is invertible and, in particular, Gjv+i.A is well-defined. Notice also 
that ([02|) yields for any Oph < tt-^/^2~^'^ '^e 

l,„n, . „o 1 
D - aphPph * -j— r 



When iV > 5, we argue exactly as in Steps [T] and [2] of the proof of Proposition 
and deduce that there exists a constant C{N), depending only on N, such that 



(5.13) WGN+iAh^nc ^CiN)\\pp4^+\ 

When iV ^ 4, our argument is different. We expand Gn+i.a as before, writting 

Gat+i^a = — aph ^ -F2j + l,A(Pph, • • • , Pph) + G6,A. 

Ar+1^2j"+1^5 

In view of dSH]) and ([ET^ (for iV = 5), this leads to 

WGN+i.Ah^nc ^ Ca^i, I M\c^' + llPphllc 

yAr+1^2j + 1^5 

In both cases, we obtain 

l|GAr+i,A||L^nc C{N) max {||pph||^+\ ||pph||S} , 
for any ctph ^ 1, so that, by (|5.8p . 

f II 1^^"^ II II® 1 

(5.14) ||GAr+i,A||L^nc G(7V)max<^ ^^i^^,^ I s: C{N,Mn 
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Finally, we consider the terms a^^™' ^^2n+i,A(i^pi , • • • ,'-^P2„+i) of sum in the 
right-hand side of (j5.9p . By (I3.7p . we have 



2n+l 



^ph 



When ^ iV, we deduce from (|4.ip that there exists a constant C(iV, A/at) such 
that \\ujp- lie ^ C{N,AfN)- On the other hand, when pj = N + 1 for some j, we can 
bound one of the norms \\ujp. \\c by ||rAr||c, and the other ones by using (|4.ip . (|5.3p 
and dEHl) to get 

llc^pjlc = M-P^\\RN\\^C{N,AfN,e)\api,\-P^, 
for Qfph ^ 1. This leads to 

^Wr-^nc ^ C'(^,A/'7v,e)|aphP""'"^max{|lrAr||c,l}, 



a. 



ph 



for Qfph 1. Combining with ([CTI)) . (j^TTll and (ISTTi)) . we conclude that 

ll^-wllL^nc C{N,J\fN, e)+(^l-e + C {N , ATn , e)\aph\^^ max{||rAr||c, 1}, 

for Qfph sufhciently small. Therefore, the norm ||7'Ar||^2nc bounded independently 
of ttph for Qfph small enough, which ends the proof of Lemma |5. II □ 

We finally prove Lemma [5T^ 

Proof of Lemma \5.Sl Given any n G {0,1,2}, it follows from recursion relations 
and (EH) that 



Therefore, given any N ^ n and — n, we deduce from (|3.3p that 

11(1 Ka - >^n)\\L-nc < ^riKl +Z^)"+*'(i'A - iy)\\L-nc 

+nK^-^\\{l+Ur+P'\UA-UML^nc- 



(5.15) 



Next, we recall that ?a — i^ls(o,2A), so that, since U{2A) = Ba, 
(5.16) 

||(i+z^)^+iz^||„^^ \\ii + u)^+^4^,^^ 

11^ ^ ^A ^IlL^nc - (i + [/(2A))^+i-"-P (1 + Ba)^+i-"-p ■ 
For the second term in the right-hand side of (I5.15[) . we use (13.21) and write 

Ua-U 



\\{l+Ur+P-'{UA-UML2nC 



^ Kq max • 



(1 + [/)W+2-n-P 

{(rTBlF+^^^'^E(2A)} IK^+^^^'^'^'lli^nc 
Since (1 -t- Ba)^+^~"~p (1 + Sa)^+^ C(iV)i;(2A), we obtain 

(l+Z^)f(i^„,A-!^„) 



L2nc 



1 — n — p 1 1 ^ ■' I 



L^nc (1-FBa)^+i-"-p 
Combining with (|5.15|) and (|5.16p . we are led to 

C{N,AfN) 



L^nc ■ 



(5.17) 



(l+Uriiyn^A-'^n) 



< 



L^nc {1 + Ba)^+^-''~p' 



for TV ^ n and Os^ps^N+l-n. 
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We next turn to the case of n ^ 3. Given any N ^ n, we assume that ()5.17p 
holds for all n ^ fc — 1 and + 1 — n, and prove it by induction for n ~ k 

and Oi^pi^N + l-k. Using (1^ and (1^ . we first infer that 



(5.18) 11(1 K,A - '^fe)|L2nc ^ 11(1 WA(j^fc-i,A - J^k-l)\ 



L^nC 



+ \Ul+U)PiUA-U)iyk- 



3^2j + l^fc 2j + i 

Y: ke = k-2j-l 

1=1 



X (^2j + l,A(l'fci.A, • • • ,l^fc2i+i,A) - ^2j + l(i^fei, ■ ■ • ji^fcsj + i)) 



L2nC 



We next estimate the first term in the right-hand side of (|5.18p using p.3p and our 
assumption. This provides 



(5.19) \\{\+UYUK(vk^i,K - i^fe-OIU^nc «i||(l +Z^)P+'(j^fe-i,A - t^fc-i)| 



L2nc 



iin I 

(1 + Sa)^^"^^''"^ IlL^nC 



For the second term, we argue as in the proof of (|5.17l) . using p. 21) and (|2.1ip : 
(5.20) 



Ua-U 



< 



(1 + ;7)Jv+2-fe-p 
(1 + Ba)^+i-'=-p' 



|(l+^)-+-'=.._r|L.,, 



Finally, we turn to the terms in the sums of the right-hand side of (|5.18p . On the 
one hand, we deduce from p.4p that 

^2j + l,A -f^2j + l — -f^2j + l,A -^2j + l,A 2j + l,A 

Hence, since p^iV + 1 — /c^A^ + l — /c^, we can apply (|3.6p and (|4.ip to obtain 



(5.21) 



{l+Uf(^F2j+i^A{'yki,Ar-- ,l^fe2, + i,A) - i^2j + l(i^fci,A, ■ • • ,i^fe2, + i,A)) 



On the other hand, the multilinearity of the function i^2j+i provides 

^2j + l(l'fei,A, • • ■ ,l^fe23 + i,A) - -F2j + l(j^/ci, ■ ■ • ,'^k2, + i) 
= F2j + i{yki,K - J^fei,l^fc2,A, • • • ,Vk2j + iA) + F2j + l{l^ki,l^k2,A -Vk2,--- ,J^fe2j + i,A) 

+ • • • + -F2J + I {Vkl ,Vk2, - ■ ■ , Vk2j I f^/C2j + l ,A — '^fe2j + l ) • 
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Therefore, we infer similarly from p.6p and (|4.ip that 
(5.22) 

(l+W)P(^F2j + i(l/fci,A,--- .Vk^^^^.K) - F2j + i(l/fei,A, • • ■ ,l'fe2, + i,A)) 
2j + l 



^c{N) J2 11(1 +wr )iic n 11(1 + ^^^''fc^.Aiic n iK^ +^)^'''=^jc 

q=l i<q t>q 

C{N,AfN) ^ C{N, Mn) 



q^'i (1 + SA)^+i-'='-f ^ {1 + Ba)n+^-^'P' 



As a conclusion, we derive from ((5?T8| . ([5J9)) . ((5?20| . (lOTj) and ((02|) that 

1 1 

Since (1 + Bp,)^+^-^-P sC (1 + Sa)^+i C(7V)A1/24, this completes the proof of 
([CT7| for n = fc. 

Notice the constant C{N,Mn) deteriorates when n increases. However, this is 
not a problem since n is limited to the set {0, • • • , iV}. Estimate (15.51) then follows 
from (|5.17p . considering the case p = 0. This concludes the proof of Lemma [5^ □ 

Appendix A. Auxiliary results on the Uehling multiplier U 

A.l. Elementary properties of U . We gather in this section some important 
properties of U , which will be useful for the proof of Lemma [A.3l in the next section. 

Lemma A.l. The function U defined in p.lOp is a non-negative, non- decreasing, 
smooth function on R-|- such that 

(A.l) C/(r) ^ and U{r) ~ ^logr. 

Its derivative U' is positive on (0, cx)), and it holds 
(A.2) U'ir) ^ A and U"{r) - 

r— >oo o7rr r— foo oTTr 

Moreover, we have 

(A.3) VreM+, -|-(l + log£;(r)) 1 + C/(r)< 1 + :^ log i;(r). 

157r dTT 

Proof of Lemma \A.l[ For the convenience of the reader, let us recall the integral 
and the explicit formulas ()2.10p of U : 

2 pi 2_ 4/0 

(A-4) Uir)^-j^^-^d. 

12-5r2 V4 + r2 „ ( ^/TTr^ + r\ 

^^^^^^^^ "'^^°gl^ V4 + r^-r j- 

Most of the statements of Lemma [A. II are direct consequences of (|A.4p . As for 
(IA.3I) . we estimate, using (|A.4p . 



CHARGE RENORMALIZATION 



21 



For the lower bound, we notice similarly that 



r2 



47r(l + ^)7o V 3 J 15^(1 + f)' 
for r E R+, so that 

(A.5) VOs^rs^l, J-(l + log^(r)) sc: ^-t^ sc: l + C/(r). 

On the other hand, we can also write 



thus when r > 1 

l + U{r) ^ -llog(l + ^) +1--^ ^ _L(i + logi5(r)). 

dTT V 2 / ovr dvr ^ ' 

The lower bound in (IA.3I) then follows from (jA.sp . □ 

A useful consequence of Lemma lA.ll is the following 
Lemma A. 2. Let be the function defined on M_|_ by 

^ ^ 1 + U{r) 

There exist three positive numbers r_ , and $0 smc/i i/iat t/ie function $ is 
an increasing diffeomorphism from (0,r_) onto (0, $0), respectively a decreasing 
diffeomorphism from (T+,00) onto (0,$o), '&~^((0,$o)) = (0,r_) U (T+,00). 
Moreover, we have 

2r 1 
(A.6) <I>(r) , and $(r' 



IStt r-i-oo rlogr 

Proof of Lemma \A.S\ From Lemma [A. 11 we see that the function $ is well-defined, 
smooth on K+, and satisfies (jA.6[) . Then we compute for r ^ 0: 

U"[r){l + U{r))^U'{rf 



$'(r) 



{l + U{r)f 



By ([XT]) and ((X2|) . we thus have $'(0) = yf^p and $'(r) --,._^oo -l/{r'^\ogr). 
Since $(0) = and $(r) as r ^ cx) by (|A.ip and (|A.2p . there exist a, 6, 5 > 
such that $ is an increasing diffeomorphism from (0, a) onto (0,5), respectively a 
decreasing diffeomorphism from (&,oo) onto (0,5). The function $ is positive on 
[a, 6], so that m = min{$(t),a ^ t ^ 6} > 0. Lemma [A. 21 follows by introducing 
$0 = min{m/2,(5}, and T_ < 1"+, the two positive numbers such that $(T'_) = 

$(r+) = $0. □ 

A. 2. A useful bound involving U . We use here results from the previous section 
to derive a bound useful for the proof of Proposition [3?2j 



Lemma A. 3. There exists a universal constant K > Q such that 



(A.7) 1 + E^^ij ^ i^logn J] (1 + U{\v,\) 

for all n ^ 1, and all {vi, ■ ■ ■ ,Vn) e (R^)". 
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If we allow K to depend on n, the optimal constant in the above inequality 
satisfies Kn — > I/Stt when n — cxd, as can be seen from the proof. The factor logn 
in (|A.7[) is therefore optimal with regard to the large-n dependence. 

Proof of Lemma \A.!^ By Lemma lA.li it holds 

(A.8) s;^^ max := J„. 

n (1 + u{\vA)) n (1 + u{\v,\)) n (1 + u{t,)) 

It is clear that taking vi — ■ ■ ■ — Vn — v shows that the maximum of the left- 
hand side of (jA.Sp is actually J„. Next, we take ti = ... = t„ = t„ in (|A.8p with 
r„ = ^157r/(nlogn). Using (|A.2p . we see that J„ > (logn)/37r for n 1. We 
will show that actually it holds J„ ^ (logn)/(37r) when n — >■ oo. In the rest of the 
proof, we assume n ^ no is such that J„ > 1. 

Let us consider a maximizing sequence {{t'f\--- ,tn'^Y\p^n for the variational 
problem defining J„. If the sequence is unbounded, then by Lemma I A. 11 

l + C/('nmaxj{4''^} 
Jn < hm 



l + U {uY&yij{tf ^) 

which contradicts J„ > 1 for n ^ rip. Therefore {t'f \ ...,tn^) is bounded in 

In this case the variational problem on the right-hand side of (jA.8|) has a maximizer, 

which satisfies the equation 



(A.9) VI « < „. m) = Xinu) ' = * I > >) h= *■ 



U 




Assume now that $1 ^ $0- By Lemma [X2l we have ^ r_ and < T+, 

for all 1 ^ fc ^ n, hence n ^ T+/T-. In particular, for n > it must hold 

< $1 < $0- Note if $1 = 0, we infer from Lemma IATTI that ii = • • • = t„ = 0, 
so that J„ = 1, a contradiction. Therefore, by Lemma lA.2[ there exist exactly two 
numbers < t„ < r_ and T„ > T+ such that $(t„) = $(T„) = $i. By (|X9l) . the 
unique possible maximizer is (r„, • • • , t„), where t„ — Tn/n e (0, T_) is such that 

(A.IO) $(t„) = $(nT„). 

The corresponding value of J„ is 

(A.11) J„- l + 



(1 + f/(T„)r 



By (|A.10[) . we must have r„ — )■ as n — )• oo. Combining (IA.10[) with (|A.6p . it 
follows that $(nr„) ~ 2T„/(157r) 0. By Lemma [A . 2 1 and since nr„ — r„ ^ T+, it 
holds riTn —7- oo. Using (jA.6[) again, we deduce that nT„ log(nr„) '--^n^oo IStt/ (2r„), 
hence finally t„ ^ -\/l57r/(ri logn). Inserting in (jA.l[) and (|A.11[) . we finally arrive 
at Jn ^ (logrt)/(37r). This ends the proof of Lemma [A. 31 □ 
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Appendix B. Proof of Proposition 13.11 

We start by showing the following lemma which provides estimates on U — Ua 
on [0,2A]. 

Lemma B.l. Let A ^ 1. For ki = 258/7r, one has 



(B.l) 



V0s$rs^2A, |[/A(r) - C/(r)| < Ki 



2£;(A) 



Proof of Lemma [Bl\ Recall that (see ([^1^ and ([XT]) ) 



(B.2) UAir)-Uir) 



in 



B(A) 



1 + 2^(1-22) 



2tt 



ZA{r) 



EiA) 



-dz 



1 



B(A) 



for ^ r ^ 2A and where 

E{A) - E{A - r) 



(B.3) 



ZAir) 



nJz.ir) (1-Z2)(1 + 2^(1-^2)) 

2A-r A 



-dz. 



r E{A) + E{A -r) E{A) ' 

We will estimate all the terms of the right-hand side of (IB. 21) . The first term is 
treated as follows, for all ^ r ^ 2A: 



(B.4) 



r 

47r 



l + ^(l-z2) 



-dz 



2 



A \ 



E{A)) " 67ri;(A)2 



Using (IB.3|) and \x\ ^ E{x), we bound the second term by 



ZTT 



i?(A) 



r-ZA(r)2 



47r(i?(A) - 2^i?(A)' 

for ^ r ^ 2 A. In order to estimate the last term of the right-hand side of (|B.2p . 
we distinguish the regions ^ r ^ A/2 and A/ 2 ^ r ^ 2 A. We calculate 



E(A) 



ZAir) il~Z^)il + rl(l-z^)) 



-dz 



i=<A) dz 



ZA(r) 



1 



■log 



1 - ZA{r) 



E{A) 



On the other hand, by (jB.3l) . 



1-^a(0 

A 



1 



r(2A-r)(A + S(A)) 



s=: 1 



6r 



1 _ ^ (^(A) + E{A - r))((A - r)E{K) + AE{A - r)) ^ ^ ' E{A) ' 

as soon as ^ r ^ A/2. Hence using log(l + a;) ^ x we infer the bound 



(B.l 



vo r A/2, 



z.M {l~z^){l+^i{l-z^)) 



-dz 



4r 



For A/2 ^ r ^ 2A, we write similarly as before 

2 _ £l 

dz ^ 



j=;(A) 

Za(.) {l-Z^){l + i{l-Z^)) 

and deduce the estimate 

(B.7) VA/2 s$ r s$ 2A, - 

TT 



3r2 



A 

E(A) 



dz 



ZA{r) 



r) {l-^Y 



3r^ 



ttE{A) ■ 
S(A)(A + ^(A)) 



B(A) 



ZaM (1-z2)(1 + 2^(1_^2)) 



-c?z 



128r 
7r£;(A) ■ 
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Estimate ([RT|) follows from jEH, ([B31) . (jRe]) and (|BJ1) . together with (|R2l) . This 
ends the proof of Lemma IB. II □ 



We now use Lemma IB. II to finish the proof of Proposition 13.11 The pointwise 
convergence of C/a when A — > oo is a direct consequence of (HH. For (|3:2|) . we first 
use (|A.3|) and (IB. II) to obtain 



VO sC r 2A, 



C/aW - U{r) 



157r\™+i 



/i57r\ 



E(r) 



(1 + C/(r))™+i 
Optimizing x (i+iogi^^i))m+i on [0, 2 A] yields 

E{r) ^ ^ r ii;(2A) 



V 2 / 2E{K){l + \ogE{r))"^+^' 



< max 



(l + log£;(r))'"+i " I ' (l + log£:(2A))^ 
Since E{2x) ^ 2E{x) for any x ^ 0, we are led to 



C/aW -C/(r) 



(1 + C/(r))"+i 



157r\™+i 



/ i07r\ 



in 



log£;(2A))™+i' £;(2A) 



1' £;(2A)}' 



On the other hand, U is non-decreasing on M+, hence, using (jA.Sp we infer 



w ^ 2A, 

' (1 + C/(r))"+i 
Using (|B.8[) . we finally obtain 

(B9) ^-^"^ 

We now recall that 

(B.IO) 

so that, for A ^ 1, 



U{r) 



(1 + C/(r))"+i 



(f) 



(1 + logi;(2A))™' 



(1 + log£;(2A))™' £;(2A)/ 



Ba = - 



dz, 



Ba^ — 

671 



-^dz = A log [iJ(A)(A + E{K))] ^ ± logi?(2A). 
1 — z dvr dvr 



Combining with (jB.9[) . we finally derive (|3.2p . We end the proof of Proposition 13. II 

□ 



by noting that (|3.3p follows directly from the definition of Ua and (jB.l 
Appendix C. Proof of Proposition 13.21 
We may define ^ (/i) by duality as follows 

for any smooth function C, and where 



in.A 



1 

2^ 



n(n 



^A 



rB 



We will use, like in [T^ p. 547], the inequality 



(C.l) 



|Tr(Q^„,AC)l 



Tr, 



Qn.ACb.P)!^^^- 



The main idea is to derive a bound of the last integral in (jC.l|) in terms of the 
norms ||(1 -|-W)"™CIIl2 and ||(1 + Z-/)^'"CIIl2+c/ , which provides an estimate of the 
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form p.6|) . by duality. The proof will depend whether we estimate the inte- 
gral in the right-hand side of (jC.ip by the norm ||(1 + U)^™'C,\\^, or by the norm 
1(1 +U)~™(\\]^2- For this reason, we split it into three steps. 

Step 1. There exists a universal constant Ci such that for all ^ 
(C.2) \\{l+UY-F:,^^{^^)\\^ ^ (Ci)"(Xlogn)- ^ j + ^)™^^. ||^ ^ 

A 2 

for all n = (^1, • • • ,Ain) e C". 

We estimate ^({p,p) as follows: 

(C.3) 

\QlZaP,p)\ ^ --irr r [ - I (/i"^b)|^(p-Pi)|/i"Hpi))x 



{2n 

n-l 

X 



X {f^^"^'\pn)\aPn-p)\ftHp))dp^ ---dp^ 



where (pj = /ij * | • and for any /3 > 0, /^'^ = tt^'VI??^ + with tt^'^ 1, if 

e ^ —1, and 7r^^^^(-) — 1|.|>a. Applying the following coroUary of (|A.7I) 

ri-l 

(1 + u{p - pn)r ^ {Kiognni + u{p ~ pi)ri[{i + u{p, - p,+^)r 

to (|C.3|) . we are led to estimating 

/ |OCCb,P)MP ^ — (iflogn)™ / Tr( n (/^^HV)V':,-Wx 

X /^+^ VzV)) (-zV)e(a;)/^) (-zV))) d77, 



where V'j = (1 + C^)™l^l for 1 ^ j s$ n, and C = (1 + C^)""ICI- Since n -I- 1 ^ 6, we 
deduce from Holder's inequality in Schatten spaces t25j, and the fact that || • He, ^ 
II • II e^, as soon as 1 ^ r ^ g ^ oo, that 



[XX 



r 1 / _ _ II 

(C.4) \Q^aP,P)\dp^Tr^K\ognr [X{W\-^^)u^ 

xf[''^'\~iV) f{'"^'\-^V)ax)f[''\-iV) 

4 Se 4 'I ©6 ^ 

We now use the Kato-Seiler-Simon inequality (see f22^ and [23] Thm 4.1]), 

(C.5) Vp ^ 2, ||/(-«V)5(a;)||6p =^ -^-xII3||lp(r3)||/||lp(r3), 

{2tt) p 
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to bound all the terms of the product in the right-hand side of (jC.4p . This provides 

f['\~tV)\h{x)\f[''\-tV) ^ #HV)|/i(x)|^ \hix)\if{''\~tV) 

4 4 Se 4 ©12 4 



■(27r) 



r 1^1 



L6 











4 


L12 


4 


L12 



for any e and e' in {—1,0,1}, and any h € L^(K^). In particular, by the critical 
Sobolev inequality, we obtain for any function h in _ff"'^(]R"^), 



(C.6) 



f^^\-^V)\h{x)\ff\~^V) ^A\\Vh\\^. 









4 




4 



L12 



for some universal constant A. Given any q ^ 2 and /? > 6/q, we then check that 



(C.7) 



LI 



du 



E{u 



= ( ) 



du 

urn WF' 



for £7^—1, while similarly. 



(C.8) 



//3 



The definition of the functions tpj gives 

(C.9) ||(l+WrVV,||i. -47r||(l+Z^r>,||c, 

which, combined with jOll, ([CJ|) . ([CJ| and JOfit, leads to 

n. 71 

/ icrc(p.p)Mp ^ logn)™ 11(1 +z^)-™ciic' n 1(1 +^r>.iic 



r 1 1 1". d-q 



for some universal constant A. When rtj — 0, we have £^(77) ("+i)/2 gfj^ ^ 
E{t])~^ drj, whereas, for = 1, 



1 1 



1 dy 



drj 



Inequality then follows with Ci = A/2 + A Eirj^^^^drj. 

Step 2. There exists a universal constant C2 such that, for n = 3 or 

{C2T{K\ogn) 



(c.io) ||(i+z^rF;A,^,„^,(j^3) 

for all fi = (^1, • • • ,^in)& C". 



nii(i+z^)">,iic> 
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The proof is similar to the proof of (|C.2|) . Since n/6 + 1/2 ^ 1, we can now 
estimate (5,\ aC(p,p) by 



(Klogn) 
2^ 



m poo 



i=2 



/h)(-zV)V„(a:)/h+^)(-zV)|| ||/^+^^(-^V)^(a;)/^)(-^V) 

4 8 Il66 II 8 2 62 / 

where the functions fjj^\ ipj and ^ are defined as in Step[TJ Using Holder's inequality 
and (|C.5|) . we can bound each norm in the right-hand side of (jC.ll[) similarly to 
(jC.6p . This provides, for instance, 



and 



f^^^\-^V)\hix)\f'f'\-^V) 



/(^"+i)(_,V)|/i(a;)|/fH-zV) 







^(.2) 








toe 




s 


LIS 


4 


L9 



T ll"llL2 



©2 (27r 

Combining with (|C.7p . (|C.8I) and (|C.9|) . we obtain 



(<^i) 



1,2 



/ \Qi^ap,p)\dp ^ A-+\Kiognr 11(1 +w)-™ciL2 n \\{i+uy 



A'j lie 



for a universal constant A. Since 



drj 



'0 E{ri) 6 +7 
for rif = 0, whereas, for — 1, 

1 1 -1 dr? 



^(77)^ 



lin I r : — r \ " 



drj 6 



.S(?7)T ATiE{r,)^ Ay Jo E{ri)i A 
we obtain (ICTOl) with C2 = + A E{r])-y^dij. 
Step 3. Lef n = 3. There exists a universal constant C3 such that 

{CsfiKlogir ' 



(C.12) 



A— 



niKi+z^r^.-ic 



/or aZZ /i = (^i,^2,A^3) € C^. 

The proof of (|C.12|) follows ideas from Ti', Section 4.3.4]. Contrarily to Steps [T] 
and[2J it relies on an explicit computation of Qg j^CiPjP) by means of the residuum 
formula for the integral with respect to the variable 77. Indeed, it holds 



(C.13) 



Q%,aC{p,p) = XI '33,aC(p,p)- 
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Here, the quantity Q'^\CiP:P) vanishes if 5 = ±(1,1,1,1), whereas, when 5 
(1,-1,— 1,-1), it refers to the expression 



Q3,aC(p,p) = 



1 



(27r)6 



^ <MP1--P2) (£3). .pO/ x^3b2-P3) (e4). ^ pO / \^ J J 

E{p) + Eip,) ''^ (P^)^-(P^) j,(p)+j;(p3) ^A b3)P-(p3)Cb3-p)dpidp2dp3, 

where -P±(p) = (^'(p) ± {ot.p + (3))/2E{p). The expression of Ql'^^C{p,p) is similar 
when (5 contains exactly one Si = I, respectively exactly one Si = —1. On the other 

hand, for S = (1, 1, —1, —1), the function Q'^^\C{PtP) is given by 



(27r) 



4''>{p)Pl{pW,ip-p^)\4''>ip,)P°{p^)x 

X l^bl -p2)|7r^^(p2)P°(p2)|^(p2 -p3)k^^b3)P°(P3)C(P3 -P)X 
/ 1 



(Eip) + Eip2)){E{pi) + E{p2)){E{pi) + E{p3)) 

! \ 

(E{p) + E{p2)){E{p) + E{p3)){E{pi) + E{p3)) 



dpi dp2 dps. 



We next estimate Q^\C{PtP) as above. For instance, when S = (1, —1, —1, — 1), 
since E{p + q) ^ E{p) + E{q) for any (p, q) E M"^, wc can compute 

|pO(p)i^(p-pi)pO(pi)| 



7:r., - 1 



3.ACiP,P)\ ^ 



dpidp2dp3 



E{p + pi)z 



^ ie2)( -P2)\ {es)f ^ l^(P2 -P3)| (^4). ^ \C{P3 - P)\ (£1). s 



£:(pi)6i;(p2)2 

so that, using (jA.7|) as above, 
(C.14) 

5lC(p,p)|dp«C (i^log3)™||A.Vl| 



E{p2)-E{p3)2 



E{p3)2E(pY 



VI II e2 



^(-,V)^3(x)^(-.V) 

i/ 2 2 



^(^2) (.3) 

^(-zV)^2(x)^(-zV) 
^(-.V)C(x)^(-zV) 
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where 



^UAP.'l)-^§^\Kip)P^-m■ 

E{p + q) 3 

The operator M^^ was estimated in Lemma 14 of 14! by ||M;pJ|g^ ^ A||V(/3i||i2, 
where A is some universal constant. By the definition of (/?!, we obtain 



(C.15) 
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As for the other terms in the right-hand side of (jC.14|) , we argue as before and, by 
(|C.14P and (|C.15p . we obtain 

(C.16) / \QSaP,p)\dp ^ ^^^!°f ^"' ||(l+^)-"CI|r n 11(1 lie, 

where A denotes some universal constant. All the terms in the right-hand side of 
(jC.13|) are similar to the one corresponding to 5 = (1,-1,— 1,-1). In particular, 
(|C.16P holds when Ql\C is replaced by Q^.aC- By duality, this completes the proofs 
of ([Cl2|) and of Step El This ends the proof of Proposition |321 □ 
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